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Fast quantum non-demolition readout from longitudinal qubit-oscillator interaction 
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We show how to realize high-fidelity quantum non-demolition qubit readout using longitudinal 
qubit-oscillator interaction. This is realized by modulating the longitudinal coupling at the cavity 
frequency. The qubit-oscillator interaction then acts as a qubit-state dependent drive on the cavity, 
a situation that is fundamentally different from the standard dispersive case. Single-mode squeezing 
can be exploited to exponentially increase the signal-to-noise ratio of this readout protocol. We 
present an implementation of this idea in circuit quantum electrodynamics and a possible multi¬ 
qubit architecture. 

PACS numbers: 42.50.Dv, 03.67.-a, 03.65.Ta, 42.50.Lc 


Introduction - Measurement in quantum information 
processors is generally realized by entangling a qubit to 
an ancillary system, with the states of the latter strongly 
depending on the qubit states. The ability to resolve 
these ancillary, or pointer, states can then result in a 
qubit measurement. This readout should be fast, quan¬ 
tum non-demolition (QND) and of high-fidelity. A com¬ 
mon approach to realize this is the dispersive regime of 
cavity QED where the electric dipole moment of an atom 
strongly couples to the electric field of a high-Q cavity [1] . 
In the dispersive regime where the qubit-cavity detun¬ 
ing A is large with respect to the coupling strength g^, 
the cavity frequency Wr is modified to take a qubit-state 
dependent value Wr ± Xj with x = 5 x/^ the dispersive 
qubit-cavity interaction. Starting in the vacuum state, 
a drive then displaces the cavity to qubit-state depen¬ 
dent coherent states |ao,i)- Resolving these two pointer 
states by homodyne detection of the transmitted or re¬ 
flected signal completes the qubit readout. Combined 
with recent advances in near quantum-limited amplifica¬ 
tion [SHU, this approach has led to high-fidelity readout 
with superconducting qubits [5HH]. 

Despite these successes, current state-of-the-art disper¬ 
sive readout suffers from several problems. First, read¬ 
out must be made faster to meet the stringent require¬ 
ments of fault-tolerant quantum computation [9]. This 
is however challenging with the dispersive qubit-cavity 
interaction taking the form iJdisp = xo^h^z- Indeed as 
illustrated by the dashed lines in Fig. [^a) under this in¬ 
teraction a coherent drive at the input of the cavity first 
pushes the two pointer states in the same direction in 
phase space before pulling them apart. There is there¬ 
fore little information about the qubit state at small mea¬ 
surement times. Second, i/disp derives in second-order 
perturbation theory from the electric-dipole interaction 
+ d)ax [I]. Since the latter does not commute with 
the measured qubit observable, Uz, dispersive readout is 
only QND in a perturbative sense. This non-QNDness 
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FIG. 1. (a) Evolution in phase space of the cavity out¬ 

put field flout for longitudinal (full lines) and dispersive cou¬ 
pling (dashed lines, x = ^/2)- Blue and red refer to qubit 
states. For the chosen parameters, the field starts at (0,0) 
to reach (0, ±1) in steady-state; the circles illustrate the posi¬ 
tion of the pointer states at characteristic times. Longitudinal 
coupling modulation and coherent drive are switched on as 
(2/7r) arctan(lOKt). (b) Pointer state separation as a function 
of time. Vertical lines correspond to the circles of panel (a). 


manifests itself with Purcell decay 7 ^ = {gx/^)‘^K [lOj . 
where k is the cavity damping rate, and with the ex¬ 
perimentally observed measurement-induced qubit tran¬ 
sitions mm- For this reason, the cavity damping rate 
cannot be made arbitrarily large and the measurement 
photon number h is typically kept well below the critical 
photon number ncrit = (^/25x)^ |13j . Increasing k and 
h would otherwise lead to faster qubit measurement. 

In this letter, we study an alternative approach that 
addresses both the slow measurement time and the non- 
QNDness. This proposal is based on a longitudinal qubit- 
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cavity interaction of the form gz{a^ + a)az- As already 
noted in Refs. [HI [15], this coupling is purely QND and 
therefore avoids Purcell decay. Here, we show that un¬ 
der appropriate driving, longitudinal interaction leads to 
an optimal separation of the two pointer states in phase 
space, without the initial slow separation that is charac¬ 
teristic of the dispersive interaction. In contrast to the 
dispersive case [isKni, we moreover find that the signal- 
to-noise ratio (SNR) of qubit readout can be exponen¬ 
tially improved by injecting a single-mode squeezed state 
in the cavity. As a possible realization of this idea, we 
discuss a circuit QED implementation based on a trans- 
mon qubit [18] strongly coupled to the flux degree of 
freedom of an oscillator. The effect of imperfections and 
a possible multi-qubit architecture are also presented. 

Longitudinal readout - Under longitudinal coupling, 
the qubit-cavity Hamiltonian reads {H = 1) 

H = uJrd'^d++ gz^zid'^ + d), (1) 

where ui^ and Wa are respectively the cavity and qubit fre¬ 
quencies, while gz is the longitudinal coupling strength. 
The realization of multi-qubit gates based on this interac¬ 
tion has already been discussed in the context of trapped 
ions [TM^ and superconducting qubits [HOES]. In 
particular, Ref. [H] proposes to modulate the bias of a 
flux qubit to realize two-qubit gates. In the absence of 
external perturbations however, this interaction leads in 
steady-state to a qubit-state dependent displacement of 
the cavity field of amplitude zLgz/{u!r + iKf2). In other 
words, longitudinal interaction is of no consequences for 
the typical case where Wj- ^ 5 zj s;. 

Here we propose to render this interaction resonant 
for readout by modulating the coupling at the resonator 
frequency: gz(t) = gz + gzCOs{uj,:t). In the interaction 
picture and neglecting fast-oscillating terms we obtain 

H = ^gzdzid'^ + d). (2) 

This now leads to a large qubit-state dependent displace¬ 
ment ± 5 z/k. Even with a conservative modulation am¬ 
plitude gz ~ 10k, the steady-state displacement corre¬ 
sponds to 100 photons and the two qubit states are eas¬ 
ily distinguishable by homodyne detection. We note that 
with this longitudinal coupling there is no concept of crit¬ 
ical photon number and large intra-cavity photon popu¬ 
lation is therefore not expected to perturb the qubit. 

While large pointer state separation can be obtained 
under a strong measurement tone in the dispersive case, 
here crucially the separation of the two pointer states oc¬ 
curs much faster. This is illustrated in Fig. [^a) which 
shows the path in phase space for both measurement 
protocols (full lines: longitudinal; dashed lines: disper¬ 
sive). The circles on this figure show the position of the 
pointer states at characteristic times until steady-state is 
reached. Clearly, with the proposed scheme the pointer 
states take the optimal path in phase space towards their 



FIG. 2. (a) SNR in units of g^l k as a function of integration 

time r. Longitudinal coupling (green line) is compared to 
dispersive coupling (dashed black line, x = '^/S) for the same 
pointer state separation, \gz\ = |e|. The dotted cyan line ac¬ 
counts for a residual dispersive coupling of x = k/ 5. The full 
brown line shows the exponential improvement obtained for a 
single-mode squeezed input state with = 20 dB. (b) Inte¬ 
gration time T required to achieve a fidelity F = 99.99 % ver¬ 
sus longitudinal coupling, (c) Cavity damping rate to reach a 
fidelity of 99.99% in 50 ns versus intra-cavity photon number 
n = {gzInY = ‘lltjKf’. Squeezing (full brown line) helps in 
further reducing the required photon number or cavity decay 
rate. The squeeze strength is optimized for each k, with a 
maximum set to 20 dB reached close to «;/27r = 1 MHz. 


maximal separation reached at steady-state. As shown in 
Fig.^b), this leads to a larger pointer state separation 
at short times. 

The consequence of this observation on qubit measure¬ 
ment can be quantified with the signal-to-noise ratio. 
The SNR is evaluated from M(t) = y/n /J” dt[aQut(t) -b 
«out( 0 ]i the measurement operator for homodyne detec¬ 
tion of the output signal dout with a measurement time 
r. From this expression, the signal is defined as \ {M)i — 
(M)o|, where the label { 0 , 1 } indicates the qubit state, 
while the imprecision noise is [(M^^(t)) -|- (./U^g(r))]^/^ 
with Mm = M — (M) [T^. Combining these two expres¬ 
sions, the SNR for the longitudinal case then reads [Mj 


SNR, = 

{--( 

1 - 

K 

tiT \ 

) 


( 3 ) 


This is to be contrasted to the SNR obtained for disper¬ 
sive qubit readout under a coherent drive of amplitude e 
and optimal dispersive coupling x = k /2 [TfllMlES] 


SNR^ = 


[i-A( 

1 —e 2 '^'^ cos iAtr) 

KT \ 

^ ) 


( 4 ) 


Both expressions have a similar structure, making very 
clear the similar role of gz and e, except for the cosine 
that is a signature of the complex path in phase space of 
the dispersive case. Importantly, for short measurement 
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times KT <C 1, we find a favorable scaling for longitudinal 
readout with SNR^, oc SNR^/kt. This advantage is illus¬ 
trated in Fig. [^a) that shows the SNR versus integra¬ 
tion time for longitudinal (full green line) and dispersive 
(dashed black line) coupling. At equivalent steady-state 
separation {gz = e), this leads to shorter measurement 
time for longitudinal coupling. This is made clear in 
Fig. [^b) presenting the measurement time required to 
reach a fidelity of 99.99% as a function of the modula¬ 
tion amplitude. 

As evidenced by the above discussion, for gz = e the ad¬ 
vantage over dispersive readout is found at short integra¬ 
tion times. This is especially true when considering the 
non-perturbative effects that affect the QNDness of dis¬ 
persive readout. As an illustration of this, Fig. He) shows 
the cavity damping rate vs photon number required to 
reach a hdelity of 99.99% in the short measurement time 
T = 50 ns. The full green line again represents longi¬ 
tudinal readout and the dashed black line the dispersive 
case. The horizontal dotted line at K/27r = 4 MHz is a 
typical value for circuit QED experiments and in partic¬ 
ular corresponds to Ref. [S] where a fidelity of 99.8% was 
achieved in 140 ns. With this k, reaching a fidelity of 
99.99% in 50 ns with dispersive readout would require as 
many as 500 photons (see vertical dotted line). For most 
circuit QED experiments, this is well above Ucrit where 
non-perturbative effects are expected to reduce the read¬ 
out fidelity. On the other hand, reaching the same goal 
with longitudinal readout only requires ~ 40 photons. 
Importantly, working at larger photon number is also a 
possibility here since there is no critical photon num¬ 
ber. Alternatively, working with ^ 10 photons requires 
a large cavity damping rate k/27t ^ 10 MHz to reach 
the above goal. Under transverse coupling, this would 
lead to significant Purcell decay. Indeed, for the typical 
value 5 x/A ~ 1 /10, the cavity-induced relaxation time 
1 /7 k ~ 1-6 gs is much smaller than current qubit relax¬ 
ation times. As already mentioned above, longitudinal 
coupling does not lead to Purcell decay [HI [15] . 

In short, the proposed approach allows to reach large 
readout fidelities in short measurement times. Reach¬ 
ing the same goal with dispersive readout requires either 
large k or large h, something that in practice would lead 
to a reduction of the readout fidelity. It is also interesting 
to point out that longitudinal readout saturates the in¬ 
equality r,^i„ > Pmeas linking the measurement-induced 
dephasing rate F^^m to the measurement rate Fmeas and 
is therefore quantum limited [24j . 

Single-mode squeezing - The SNR of longitudinal read¬ 
out can also be exponentially improved with a single¬ 
mode squeezed input state on the cavity. For this it suf¬ 
fices to chose the squeeze axis to be orthogonal to the 
qubit-state dependent displacement generated by gzit). 
In Fig. Ha), this corresponds to orienting the squeeze 
axis along the vertical axis. With this choice, and since 
the squeeze angle is unchanged under evolution with lon- 





FIG. 3. (a) Circuit QED implementation of longitudinal cou¬ 
pling with a transmon qubit of Josephson energies Eji = 
Ej{l d)/2, Ej 2 = Ej{l — d)/2 with d £ [0,1]. (b) Qz and 
gx versus flux. Around <l>x = 0, Qz as a linear dependence 
with flux. The spurious transverse coupling px results from 
qubit asymmetry, here d = 0.02. The dashed lines represent 
the asymptotic expressions Eqs. (§ and Q, the full lines to 
exact numerical result [24|. (c) Transmon frequency versus 
flux for Ej/h = 20 GHz, Ej/Ec = 67 and d — 0.02. 


gitudinal coupling, the imprecision noise is exponentially 
reduced and the signal-to-noise ratio simply becomes 
e’'SNRz, with r the squeeze parameter [ 21 ]. This ex¬ 
ponential enhancement is apparent from the full brown 
line in Fig. [^a) and in the corresponding reduction of 
the measurement time in Fig. Hb). 

This is in stark contrast to standard dispersive read¬ 
out where single-mode squeezing can lead to an increase 
of the measurement time [HidT]. Indeed, under disper¬ 
sive coupling, the squeeze angle undergoes a qubit-state 
dependent rotation. As a result, both the squeezed and 
the anti-squeezed quadrature contributes to the impreci¬ 
sion noise. We note that the situation can be different 
in the presence of two-mode squeezing m where the 
present exponential increase in SNR can be recovered by 
engineering the dispersive coupling of the qubit to two 
cavities m- 

Circuit QED implementation - We now turn to a pos¬ 
sible realization of this protocol in circuit QED. Longi¬ 
tudinal coupling of a flux or a transmon qubit to a LC 
oscillator was already discussed in Refs. [IHCl]- There, 
longitudinal coupling results from the mutual inductive 
coupling between a flux-tunable qubit and the oscillator. 
As another example, we follow the general approach de¬ 
veloped in Ref. [26] and focus on a transmon qubit that 
is phase-biased by the oscillator. Fig.^a) schematically 
represents a lumped version of this circuit. In practice, 
the inductors can be replaced by a junction array 
both to increase the coupling and to reduce the size of the 
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qubit’s flux-bias loop. An in-depth analysis of an alter¬ 
native realization based on a transmission-line resonator 
can be found in Ref. [24] . 

The Hamiltonian of the circuit of Fig.|^a) is similar to 
that of a flux-tunable transmon, but where the external 
flux <i)x is replaced by <l>x -I- <5 with S the phase drop at 
the oscillator |2S]. Taking the junction capacitances to 
be equal and assuming for simplicity that 1 

with Zq = yjLjC and Rk the resistance quantum, this 
Hamiltonian can be expressed as H = + Hq + Hqr- In 

this expression, = uJrO'a is the oscillator Hamiltonian 
and Hq = is the Hamiltonian of a flux-tunable 

transmon that we write here in its two-level approxima¬ 
tion [THj. The qubit-oscillator interaction takes the form 
tiqr = + a)<5’a; + <7z(a^ + a)o'z with [21] 



and where Ej is the mean Josephson energy, d the 
Josephson energy asymmetry and Eq the qubit’s charg¬ 
ing energy. Expressions for these quantities in terms of 
the elementary circuit parameters are given in Ref. |24j . 
As desired, the transverse coupling ^x vanishes exactly for 
a symmetric transmon with d = 0, leaving only longitudi¬ 
nal coupling g^. Thanks to the phase bias rather than in¬ 
ductive coupling, gz can be made large |5^. For example, 
with the realistic values Ej/h = 20 GHz, Ej/Ec = 67 
and Zq = 50 fl we find gz/27r « 135 MHz x sin (tt^x/^o). 
The flux dependence of both g^ (blue lines) and g^ (red 
lines) for a finite asymmetry d = 0.02 are illustrated in 
Fig-ib). The dashed lines refer to the above asymptotic 
expressions for g^ and g^ while the full lines are exact nu¬ 
merical results. Modulating the flux by 0.05<i>o around 
$x = 0, we find ^z/Stt ^ 21 MHz. This is accompanied 
by a small change of the qubit frequency of ^ 40 MHz, 
see Fig. |^c). Importantly, this frequency change does 
not affect the SNR under longitudinal readout [2i] . 

Tolerance to imperfections - When d ^ 0, a finite 
transversal coupling is present. This is illustrated in 
Fig. §b) where for a realistic value of d = 0.02 [21] 
and the above parameters we find (7x/27r « 13 MHz x 
cos (7r$x/4>o). The effect of this unwanted coupling 
can be mitigated by working at large qubit-resonator 
detuning A where the resulting dispersive interaction 
X = dx/^ can be made small. For exemple, the above 
numbers correspond to a detuning of A/27r = 3 GHz 
where 5.6 kHz. It is important to emphasize 

that, contrary to dispersive readout, the longitudinal ap¬ 
proach is not negatively affected by a large detuning. The 
large detuning moreover reduces Purcell decay which, for 
the above g^, can very easily be made small enough to 
avoid the qubit from being Purcell limited. 


When considering higher-order terms in Zq/Rk, the 
Hamiltonian of the circuit of Fig. [^a) contains a 
dispersive-like interaction Xz^^dd^ even at d = 0. For 
the parameters already used above, we find Xz/27r ~ 5.3 
MHz [ 23 ]. Contrary to the standard dispersive coupling, 
Xz cannot be made small by detuning the qubit from 
the resonator. However since it is not derived from a 
transverse coupling, it is not linked to any Purcell de¬ 
cay. Moreover, at small integration times SNR^ is not 
affected by a finite dispersive-like coupling [ 23 ]. This is 
illustrated in Fig. [^a) where the cyan dotted line, cor¬ 
responding SNRz in the presence of x = k/ 5, is barely 
distinguishable from the ideal case. 

Finally, when discussing the exponential gain in SNR 
provided by single-mode squeezing, we have assumed a 
source of broadband pure squeezing. The effect of a finite 
squeezing bandwidth F was already studied in Ref. m 
and only leads to a small reduction of the SNR for F 3> k. 
On the other hand, deviation from unity of the squeezing 
purity P leads to a reduction of the SNR by 1/ VP- The 
SNR being decoupled from the anti-squeezed quadrature, 
the purity simply renormalizes the squeeze parameter. 

Multi-qubit architecture - A possible mutli-qubit ar¬ 
chitecture consists of qubits longitudinally coupled to 
a readout resonator (of annihilation operator Oz) and 
transversally coupled to a high-Q bus resonator (ox). The 
Hamiltonian describing this system is 

H = UJrzdldz Wrxoj^Ox + ^ 

^ (7) 

+ 5l5zjxtz,(at+az) + ^ ^x)- 

j 3 

In this architecture, readout is realized by taking advan¬ 
tage of the longitudinal coupling while logical operation 
are realized using the bus resonator. A scalable archi¬ 
tecture taking advantage of longitudinal coupling is dis¬ 
cussed at length in Ref. m- Here, we again consider 
that the longitudinal coupling of each qubit can be mod¬ 
ulated independently. We take this modulation to be 
9zj(t) = ffz + 5z cos(a;rt -I- (fij) where the phase ipj is ad¬ 
justable. In the interaction picture and neglecting fast- 
oscillating terms, the longitudinal coupling becomes 

-ffz = -bh.c. (8) 

3 

This effective drive on the resonator displaces the field 
to a multi-qubit-state dependent coherent state allowing 
single-shot multi-qubit measurements. For two qubits, 
tpj = j7r/2 leads to four well separated states in phase 
space. Other choices of phase lead to overlapping pointer 
states corresponding to different multi-qubit states. Ex- 
emples are ipj = 0 for which the two-qubit states |01) and 
110) are indistinguishable, and ipj = jir where these states 
are replaced by |00) and |I1). This can be exploited to 
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create entanglement by measurement [30) . The 3-qubit 
GHZ state is obtained with ipj = j27r/3 [24]. 

Conclusion- We have shown that modulating longitu¬ 
dinal coupling between a qubit and an oscillator leads to 
fast QND qubit readout. Because of the optimal motion 
of the pointer states in phase, the measurement time is 
reduced with respect to standard dispersive readout. For 
the same reason, this approach can be further improved 
by using single-mode squeezing. 
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I. SIGNAL-TO-NOISE RATIO 

We derive the signal-to-noise ratio (SNR) for different qubit readout protocols: Longitudinal coupling in Sec. I A; 
Longitudinal coupling with single-mode squeezed states in Sec. IB; Standard dispersive coupling in Sec. IC; Longi¬ 
tudinal and transverse couplings in the dispersive regime in Sec. ID. 


A. Longitudinal coupling 

1. Modulation at the resonator frequency 

We consider a qubit longitudinally coupled to a resonator with the Hamiltonian 

H = oOrO^a-i- iwadz -b [9z{t)a^ -P g*{t)d]az. (SI) 

In this expression, is the resonator frequency, Wa the qubit frequency and the longitudinal coupling that is 
modulated at the resonator frequency: 


gz{t) = + Iffzl cos(wrt -b ip). 


(S2) 


In the interaction picture and using the rotating-wave approximation (RWA), the above Hamiltonian simplifies to 

H = lig^ci'^ -\-g*d]d^, (S3) 


where gz = Iffzle*''’ is the modulation amplitude. From Eq. (S3), it is clear that the modulated longitudinal coupling 
plays the role of a qubit-state dependent drive. The Langevin equation of the cavity field simply reads 

d =—i^gzCTz — \nd — y/Hdi-a, (S4) 

where din is the input field [SI]. Taking this input to be the vacuum, the input correlations are then defined by 
(din(t)d;jj(F)) = [din(t), d|jj(f')] = 5(t — t'). Using the input-output boundary condition dout = din + '/nd [SI], 
integration of the Langevin equation leads to 


O^out (t) 





Pt ^ ^ 

dont{t) = din{t) -K ^din(t'). 

J —oo 


(S5) 


where Oout = (dout) stands for the output field mean value and dout = dout — ctout its fluctuations. Because here the 
qubit-dependent drive comes from modulations of the coupling, and not from an external coherent drive, there is no 
interference between the outgoing and the input fields. As a result, a = ctout/v^ and the intracavity photon number 
evolves as 


{d)d) 



1 

I - e 2 



(S6) 


The measurement operator corresponding to homodyne detection of the output signal with an integration time r 
and homo dyne angle (j)\, is 


M(r) = yK / dt[dj)„t(t)e*‘^*'-bdout(i)e *‘^‘']. 
Jo 


(S7) 
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The signal for such a measurement is (M) while the noise operator is = {M). In the presence of a qubit, the 
measurement signal is then 


(M)i 


(M)o = 4|g,,|sin((^- (^h)'r 


r 2 / 

1 \ 

KT \ 

1 _ e-5'^^j 


(S8) 


On the other hand, the measurement noise is equal to (M^(t)) = kt. Combining these two expressions, the signal- 
to-noise ratio (SNR) then reads 


SNR2 




8|5z 


■ sin^(i^ — fph) KT 


1 - 

nr \ 


2 


(S9) 

(SIO) 


The SNR is optimized by choosing the modulation phase (p and the homodyne angle such that (p — (ph 
With this choice, the optimized SNR finally reads 


SNR = 

n 



^ mod TT. 


(Sll) 


At long measurement times (r ^ ^/k), the signal-to-noise ratio evolves as SNR = while in the more 

experimentally interesting case of short measurement times we find SNR = In short, the SNR increases 

as much faster than in the dispersive regime where the SNR rather increases as (cf. Sec. IC). 


2. Measurement and dephasing rates 

Following Ref. [S2], to evaluate the measurement-induced dephasing rate we apply a polaron-type transformation 
on Hamiltonian Eq. (S3) consisting of a displacement of d by —ig^dzlK. Under this transformation, the cavity decay 
Lindbladian K'D[d\f} = 'D[d]p = apa^ — ^{a^a,p}, leads to z\p where r,^i„ = 2|pzp/«; is the measurement- 

induced dephasing. On the other hand, the measurement rate is obtained from the SNR as Tmeas = SNR^/(4 t) = 
2|pzp/«; [S3]. The relation between the dephasing and the measurement rate is then Fmeas = This is the bound 
reached for a quantum limited measurement [S3]. 


3. Modulation bandwidth 


We now consider the situation where the longitudinal coupling is modulated at a frequency ujm ^ tOy, i.e. gz(t) = 
5z + |Pz| cos(uJmt + p). Assuming that the detuning = tOm — Wr is small with respect to the modulation amplitude 
the Hamiltonian in a frame rotating at the modulation frequency now reads under the RWA as 

H =-Ama^d+llgzO^ + g*d]dz, (S12) 

The corresponding Langevin equation is then 

d =—i^gzdz + iAmd — ^Kd — ^/Kdin, (S13) 


yielding for the output field. 


C^out(^) — 


*5z 


K — 2iAr, 


1 - 


dout{t) = din{t) - K dCe^*'^™ 20(‘ 


(S14) 


From these expressions, the measurement signal is then 

4|5z| 


(M)i - (M)o = 


v'i + (^) 

8|gz|/K 

l + (^) 

8|gz|/K 

l + (^) 


sin[(p — 4>-h + arctan(2Am/«^)] t 


2 sin[ip — (fih + 2arctan(2Am//«)] 

2 sin[ip — <))h + 2arctan(2Am/'«) + A^r] 


(S15) 
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While the signal is changed, the noise is however not modified by the detuning. From the above expression, given a 
detuning Am and a measurement time r, there is an optimal angle tp that maximizes the SNR. 


B. Longitudinal coupling with squeezing 


We now consider the situation where the modulation detuning is zero and where the input field is in a single-mode 
squeezed vacuum. This leaves the signal unchanged, but has we now show leads to an exponential increase of the 
SNR with the squeeze parameter r. Indeed, in the frame of the resonator, the correlations of the bath fluctuations 
are now 


f (dln(i)din(i')} (din(t)djn(t')}] ^ f sinh^ r i sinh2r 

\(dl(i)dl(t')} (din(t)dl(t')}J ^sinh2re-2*« coshV J ^ 


(S16) 


where we have assumed broadband squeezing with a sqneeze angle 0. The measurement noise is then 


{Mn{t)) = k dt dt' + {dont{t)dl^tit')) 


+ {d 

out it)d 

out (t'))e 




{dl,{t)dL{t'))e^^‘^- 


The output-field correlations are easily obtained from 

[ {douMdontit')) {dout{t)doutit'))\ ^ / {dln{t)din{t')) {(knit)din{t'))\ 

\{dLtit)dLtit')) (dout(0rfLt(^'))/ \{dlnii)dlnid)) (rfin(t)4n(^'))/ 

which holds here since the drive is ‘internal’ to the cavity. As a result 

(M^(t)) = {cosh(2r) + sinh(2r) cos[2((/)h — 0)]}kt. 


(S17) 


(S18) 


(S19) 


The noise is minimized by choosing 0 according to 0 — = f mod tt. With this choice, the SNR reads 

SNR(r) = e’’ SNR(r = 0). (S20) 

The SNR is thus exponentially enhanced, leading to Heisenberg-limited scaling as discussed in Ref. [S4]. 


C. Dispersive coupling 


For completeness, we derive here the SNR for dispersive readout; this result can already be found in Ref. [S2]. In 
the dispersive regime, the qubit-cavity Hamiltonian reads 

H = uj-rd) d + ddz, (S21) 

where x = 9^1 ^ is the dispersive shift. The Langevin equation of the cavity field in the interaction picture then reads 

d =—ixdzd — \Kd — \/Kd'm. (S22) 


With a drive of amplitude e = |e|e*‘^‘* on the cavity at resonance, the input field is defined by its mean = {dm) = 
and fluctuations din = din — ce-m- Integrating the Langevin equation yields 


i^out (^) /— ^ 


^ -iVqhi^z) 


X _ 2e“(A(a^> + 5K)t-i-5i<Pqb(a^> 


dout(t) = din(t) -K f dt'e-(*>^<"^> + 5'=)(‘-‘')din(t'), 

J —oo 


(523) 

(524) 


where ipqh 
is as 


2arctan(2y/K) is the qubit-induced phase of the output field. Moreover, the intracavity photon number 


2|e| 


(d'^’d) = ( TlLi j cos^dv^qb) 1 —2cos(xt)e + e 


(S25) 
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From the above expressions, the measurement signal is 


Ms,\i) - Ms^\o) = 4|e| sin((/9qb) sin((,i)d - M t 


1 


4 

KT 


COS^(i(^qb) 


sin(xr + v?qb) 
sin((^qb) 



(S26) 


On the other hand, the measurement noise is simply equal to (-M^(r)) = kt. The measurement signal is optimized 
for (/3d — </'h = f mod tt and at long integration times by i/3qb = f) or equivalently x = [S2]. For this optimal 

choice, the SNR then reads 


SNR = 


1 - 


1 —e cos kKT 


KT 


(S27) 


SNR=:^^(«r)5/2. 


At long measurement times, the SNR evolves as SNR = \/8 — ^/kt and at short measurement times it starts as 

1 I 

vTs 


D. Effect of a residual transverse coupling 


We consider the presence of a spurious transverse coupling in addition to the longitudinal coupling g^, 

H = + Wa(T^ + {[^x + 5x cosiuj^t + Vk)]^x + [ffz + ffz cos(wrf + + a). (S28) 

We now assume that g^ ^ A and we follow the standard approach to eliminate the transverse coupling [S5]. To 
leading order in g^/A and under the RWA, we hnd in the interaction picture 


H =\{x^- 2xxz)(72 + + \ [g^a) + g*a], 

with the dispersive shifts x = Xx - 4xxz, Xx = gi/A and Xxz = ffxffz/A. 

Going to an interaction picture also with respect to the first term of Eq. (S29), our starting point is 

H = xa'^aa-z + ^Ig^.a'^ + gta]^z- 

This leads to the Langevin equation 

a = -i\g^a^ - {ix + 5^)0 - y/nain- 
Following Sec. I A, the measurement signal is 


(M)i - (M)o = 4|g,,|sin((/3-(/)h)rcos^(iv3qb)<i 1-cos(v3qb) - cos((,3qb + X'r)e 

KfT L 


(529) 

(530) 

(531) 

(532) 


where as before we note the dispersive-coupling-induced rotation (/jqb = 2arctan(2x/K). Again as above, the mea¬ 
surement noise is not changed by the dispersive shift. Choosing ip — (j)^^ = ^ mod tt, the SNR finally reads 


SNR(x) = v^—V^cos^dv^qb) <1 1 - — I cos((/?qb) - cos((/3qb -b XT)e ^ 


2 r 


KT L 


(S33) 


The residual dispersive coupling reduces the value of the SNR, with the decrease behaving differently at long and 
short measurement times. At long measurement times, the dispersive coupling reduces the SNR by 

SNR(x) cs cos^(4(/3qb)SNR(x = 0) = ^^^^^^^SNR(x = 0), for r »!/«;. (S34) 

The SNR is not affected for x ‘C k/2. Interestingly, at short measurement times the SNR is completely independent 
of the spurious dispersive shift to leading orders 

SNR(x) cs — i/ir), forr<Cl/«;. (S35) 

In short, the SNR is not affected by a spurious transverse coupling for short measurement times t <C 1/k. 
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II. CIRCUIT QED REALIZATION 

We now turn to a possible realization of longitudinal coupling in circuit QED. While a lumped circuit is presented 
in the main text, here we focus on a transmon qubit that is phase-biased by a coplanar waveguide resonator. We 
recover the lumped element results in the appropriate limit. We emphasize that the numerical results presented here 
are obtained for this coplanar realization. For this reason these numerical values differ from, but are compatible with, 
what is found in the main text. 


A. Circuit Lagrangian and Hamiltonian 


As illustrated in Fig. SI, we consider a transmon qubit coupled to the end of the center conductor of a A/4 resonator. 
To increase the coupling a Josephson junction can be inserted in the center conductor of the resonator at the location 
of the qubit; see dark region in the resonator’s center conductor of Fig. SI. This is reminiscent to the approach to 
ultra-strongly couple a flux qubit to a resonator discussed in Ref. [S6]. The modelling of the present circuit closely 
follows that reference and we include the details here for completeness. 

The Lagrangian of this circuit, £ = £„ + £q + £qr, is composed of three parts consisting of the bare resonator £r, 
qubit £q and interaction £qr Lagrangians. From standard circuit theory [S7], the resonator Lagrangian takes the form 




C° -r 


1 


C'jr ;2/ 


j dx + Ejrcos[ilj{0)] + ( 0 ). 


(S36) 


where ipi^x) is the position-dependant field amplitude inside the resonator. In this expression, we have assumed that 
the resonator has total length L with capacitance and inductance per unit length. In the single mode limit, 
we write %l){x,f) = ip{f)u(x) where u{x) is the mode envelope [S8]. The Josephson junction in the resonator’s center 
conductor has energy Ej,. and capacitance Cjr. This junction creates a discontinuity '0(0) ^ 0 in the resonator 
field that will provide the desired longitudinal interaction. The coupling inductance can be replaced by a SQUID, or 
SQUID array, without significant change to the treatment [S9]. This Lagrangian was already studied in Refs. [S6, S8]. 

The transmon qubit is composed of a large shunt capacitor Cs and of Josephson junctions of energies Eji and Ej 2 , 
and capacitances Cgi and Cg 2 respectively. In terms of the branch fluxes defined on Fig. Sib), the qubit Lagrangian 
taking into account the coupling to the resonator is 


£ +£ 


^91 A2 I ^92 12 I A2 


-k Eji cos[0i] -k Ej 2 COs[ 02 ]- 


(S37) 


With the resonator phase bias A0 across the coupling junction, the flux quantization around the qubit loop reads 
(5 = 01 -k 02 = ‘f’x + 0(0). Defining 0i = 5/2 — 6 and 02 = 5/2 -k 9, we can write the above as 


Aq -kA 


qr — 


/ Cql+^//lfi^2 

V 27r y [ 2 


Cs + Cgl + Cg 2 ^.2 ^ 


+ Eji cos[J/2 -9]+ Ej 2 cos[J/2 -k 9]. (S38) 


Assuming the resonator phase-bias to be small 0(0) <C 1, we find to zeroth-order in 0(0) the usual Lagrangian of 
an asymmetric flux biased transmon qubit [SIO] 

‘ho \ Cs + 

^ J 


Cgl + C, 


92^2 


9"‘ + Ejy, [cos(<i>x/2) cos 9 — d sin($x/2) sin 9 ], 


(S39) 




FIG. SI. a) Transmon qubit coupling longitudinally to a resonator. A Josephson junction at the location of the qubit can be 
used to increase the coupling, b) Closeup of the qubit with the dehnitions of the branch fluxes used in the text. 
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where Ejy, = Eji + Ej 2 and d = {Ej 2 — Eji)/Ej^ is the junction asymmetry. Defining n to be the conjugate charge 
to 9, the corresponding Hamiltonian is 


= AEcn^ 


Eji: 


cos($x/ 2 ) cos 6 


(isin($x/ 2 ) sin 0 


(S40) 


with the charging energy Eq = l2{Cs + Cqi + Cq 2 ). 

For clarity, we project on the qubit subspace {|0), |1)} where the total Hamiltonian takes the form 

H = w-co)d + K{a)of' -\--^dz-\-Hqj-- (S41) 

In this expression, Wa is the qubit transition frequency and K the Kerr non-linearity. The latter can be made small 
and in particular negligible with respect to the photon decay rate k. To first-order in V’(O) sad in the same two-level 
approximation, the interaction Hamiltonian reads 


where 


Hm = EjIZ 


V’(O) 


qr 2 

= 9z{a^ + d)a^ + g^ia^ + a)fo 


sin($x/2) cos 9 + d cos(<l)x/2) sin 9 


Ejy, 

ffz($x) = - V’rms^—[wil($x) “ TOOo(^x)], 
Eji 

^x(4^x) —V^rms 2 aroi(4*x)- 


(S42) 


(543) 

(544) 


To obtain these expressions we have used tp^O) = u{0)tprms{d'^ +a), with i/'rms = \/d'KZ^/RK- Here is the resonator 
mode impedance and Rk = hje^ the quantum of resistance. The mode amplitude m( 0) at the location of the 
resonator’s junction is related to the participation ratio of the coupling inductance = Ejr{2'iT /to the total 
inductance of the resonator mode as m(0) k -q = Lj/Lr [S 8 ]. The lumped element limit discussed in the main text 
is obtained for 77 —)• 1. Moreover, we have defined = {i\ sin($x/2) cos0-|-(icos(<l)x/2) sin0|j). Eqs. (S44) and (S43) 
are calculated numerically by diagonalizing the transmon Hamiltonian; see full lines in Fig. 3(b). 


B. Numerical evaluation of the coupling strength 

We now present a set of possible parameters for this circuit. Is is very important to emphasize that these numbers 
relate to the eoplanar architecture discussed here and not to the lumped-element version found in the main text. We 
consider a A/4 resonator of characteristic impedance Z^ ~ 50 O, total length 3.2 mm and capacitance per unit length 
= 0.111 nF/m and inductance per unit length if = 0.278 /rH/m. The coupling inductance consists of an array 
of 13 josephson Junctions with a total Josephson energy Ej^jh = 420 GHz. Following Ref. [S 8 ], we find the first 
resonator mode to be at frequency Wr/27r = 9.99 GHz, with a characteristic mode impedance Zr = 44.7 H. The 
participation ratio of the array is found to be 77 ~ 0.53. With these parameters, the resonator bias on the qubit is 
Ls« 0.079. 

The transmon Josephson and charging energies are Eji^/h = 20 GHz and Ec/h = 0.3 GHz respectively, and 
the asymmetry d = 0.02 [Sll]. These parameters yield a;a/27r = 6.6 GHz when evaluated at the flux sweet-spot 
$x = 0- The magnetic flux is modulated around this flux value with a maximum excursion of <l>x/27r = 0.1. At 
the sweet-spot, the qubit-resonator is A = Wa(0) — ojr ^ 2 tt x 3.4 GHz. The longitudinal coupling Tjz and the 
various spurious interactions found numerically using these numbers are summarized in Table I. In particular, we 
And 5z/27r = 21.5 MHz and a small transverse coupling resulting in a negligible dispersive shift Xx = 5x/^ ~ 13 KHz. 
The maximum change in qubit frequency is Wa(0) — a;a(‘l’x) ~ 27r x 170 MHz. 

The array of N junctions allows for strong longitudinal coupling while reducing the resonator nonlinearity. Gom- 
pared to a single junction of equal energy, the non-linear Kerr effectively decreases as K = Ki/N"^ [S12] and from the 
parameters above we get Ar/27r = 48 KHz using 13 junctions. 
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TABLE I. Longitudinal coupling rate 'Qz and spurious leading-order and second-order couplings for a transmon coupled to a 
A/4 resonator with a coupling inductance consisting of an array of 13 Josephson junctions. See text for the parameters. The 
numbers refer to the coplanar architecture, not the lumped element version discussion in the main text. 


Couplings 


Rates 

(27rx MHz) 

Longitudinal [Eq. (S46)] 

9^ 

21.5 

Transverse [Eq. (S47)[ 

flx 

6.56 

Dispersive 

5^/A 

0.013 

Non-linear dispersive [Eq. (S54)j 

2Az 

2.6 

Non-linear transverse [Eq. (S55)[ 

Ax 

0 

Kerr nonlinearity 

K 

0.048 

Qubit frequency shift (from flux-drive) [Eq. (S52)] 


4.4 X 10“® 


C. Asymptotic expression for the coupling strengths 

To obtain asymptotic expressions for the qubit-resonator coupling, we consider the transmon as a weakly anharmonic 
oscillator [SIO]. In this situation, 9 « 




2Ec 


1/4 


cos(<I>x/2) {i\ (6'l' +b)\j), rrii. 


2Ec V/"sin($x/2) ,. 




z|(St+6)2|z). (S45) 


Restricting to the {|0), |1)} subspace and using Eqs. (S44) and (S43), we hnd the asymptotic expressions 


Ejy, 


2Ec 


1 1/2 


5x ~ dEjY: 


2Ec 


TT Z v . , / 

— sm($x/2) r], 
Kk 


1/4 


TT Zr- , / ^ \ 

— cos($x/2) Z7- 
Kk 


(546) 

(547) 


These correspond to Eqs. (5) and (6) of the main text in the lumped-element limit zy —)■ 1. 


D. Upper bounds on the coupling strength 


An upper bounds for the longitudinal coupling is obtained by expressing the (/z in units of the qubit frequency. 
Using Eq. (S46), we find 


^ ^ $X / TT^r 

Wa 8 V ELk ^ 


(S48) 


Eor Zr ~ 50 O resonator, maximal coupling is reached in the lumped-element limit of the oscillator where zy —)■ 1 and 
we get (/z/wa ~ 0.006, or about 4 times the coupling we obtained in our previous example giving 5z/27r ~ 42.4 MHz. 

Even larger values of y/z can be achieved by using lumped LC-circuit comprised of a superinductance of large 
impedance ~ i?if/4 [S13]. Eor participation ratios in the range zy ~ [10“^, 1], the coupling is enhanced by a factor 
of ~ 10 to 5z/wa = 4>xZyv^/8 ~ [7 x 10“"^, 7 x 10“^]. While the previous circuit model would have to be refined to 
take the large phase bias Az/ into account, it is safe to say that the larger and more compact the coupling inductance 
is, the stronger the longitudinal coupling will be. 



III. SPURIOUS COUPLINGS AND IMPERFECTIONS 


A. Effect of the flux drive 


From the first term in Eq. (S38), an oscillating external magnetic flux $x(0 = ‘f’xCOs(a;rf + t/^) at frequency ujd 
leads to an effective voltage drive on the resonator 


C 


r,drive — 



(C'gl + C'g2)‘i’x'*/’(0). 


(S49) 


This term leads to an effective drive on the resonator Hr,d = CrS Qf drive amplitude Cr 


Cr 


I I 

V’rms 

otLcj 


(S50) 


where Ecj = e^/[2(C'qi + Cq 2 )]- With the above circuit parameters, we hnd er/27r ~ 10 MHz. If desired, the effect 
of this drive can be cancelled by an additional on the input port of the resonator. Otherwise, this simply leads to an 
additional qubit-state independent displacement of the cavity held that does not affect the signal-to-noise ratio. 

In the same way, flux modulation also directly drives the qubit. This is caused by the last charging energy term of 
Eq. (S38) and yields 


Hq,d = 



{Cq2 — Cql)^xq, 


(S51) 


From the asymptotic expression qK, i [Ejli2Ec)^^^ (&'!' — h) we get Hq^d = CqS _|_ d.c. where 


e 


q 


dCj / Ej 

U 2 EC 


1/4 


$xWr 


(S52) 


is identically zero for symmetric junctions {d = 0). With the above parameters and CjICy. = 0.01, we hnd that 
eY^j^-K ~ 0.36 MHz. Given the strong qubit-resonator detuning of several GHz, this is however of no consequences. 


B. Higher-order interaction terms 

To second order in At// in Eq. (S38) yields the additional interactions 

= (a'l' -f a)^(A,j;(Ta; -f A^dz), (S53) 


where 


((0| cc - d ss|0) - (1| cc - d ss|l)), (S54) 

Id 

(1| [cc - ss]|0), (S55) 

and where we have dehned cc = cos($x/2) cos0 and ss = sin($x/2) sin0. 

Under hux modulation, Afc(t) = A^ -|- A^ cos(wrt + V^r) with k = x,z. With the RWA this leads to 

« {2d^d + l)(K^ax + A^d^). (S56) 

For d small and hux modulations around $x = 0 then A^, = 0. On the other hand, the dispersive-like interaction of 
amplitude Xz = 2 A 2 has a magnitude of Xz/27r ~ 2.6 MHz for the above circuit parameters. In the lumped-element 
limit with Zy = 50 H, we rather hnd Xz/27r = 5.3 MHz. Since this dispersive-like interaction originates from a 
second-order correction, we always hnd (/z/Xz ~ 10- In practice, we can therefore chose Xz < « such that this spurious 
coupling will not affect the SNR at short integration times. 
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FIG. S2. (a) Multi-qubit architecture consisting of qubits longitudinally coupled to a readout resonator (dark blue) and 

transversally coupled to a bus resonator (brown), (b)-(e) Phase space representation of the pointer states for two and three 
qubits, (b) Two-qubit joint readout {ipj = j-n/2), (c-d) Two-qubit entanglement by measurement. The green circles correspond 
to the Bell states 1^'’*') = ^[|0> 1) + l^i 0)1 1“^^) = 0) + l^> 1)]’ obtained for (pj = 0 and ipj = Jtt, respectively, (e) 

The green circle corresponds to the GHZ state |GHZ) = ^[|0,0,0) -I- |1,1,1)]. This configuration is obtained for (pj = j2'Kl3. 


IV. MULTIPLE QUBIT READ-OUT IN (X-|-Z) CQED 

Fig. S2(a) illustrates a possible multi-qubit architecture: qubits (green) are longitudinally coupled to a readout 
resonator (dark blue) and transversally coupled to a bus resonator (brown). The readout resonator is used only when 
the longitudinal couplings are modulated and individual flux control for each qubit is possible with separated flux 
lines. The system Hamiltonian reads 

H = -j- CUrxdx^x T ^ ^ ^ ^ ^ 9zj^zj{^\ T dz) T ^ ^ “t” dx). (S57) 

3 3 3 

In the interaction picture and neglecting fast-oscillating terms, the longitudinal coupling becomes 

-ffz = (^l^z +h.c. (S58) 

3 

The position of the qubits and the coupling inductances can be adjusted to get equal longitudinal coupling 
strengths \gzj\- By properly choosing the phases ipj it is possible to operate in the joint qubit readout mode 
or in the entanglement by measurement mode. Examples for 2 and 3 qubits are shown in Fig. S2(b)-(e): joint 
measurement of two qubits, synthesis of Bell states and 3-qubit GHZ state by measurement, cf. caption. 
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